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Abstract 

The stability problem for the O(N) nonlinear sigma model in the 2 + e di- 
mensions is considered. We present the results of the 1/N 2 order calculations 
of the critical exponents (in the 2 < d < 4 dimensions) of the composite op- 
erators relevant for this problem. The arguments in the favor of the scenario 
with the conventional fixed point are given. 
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In a recent time the possibility of a scenario with the nonconventional fixed points in 2 + e 
expansions for various models was widely discussed. It has been firstly observed by Kravtsov, 
Lerner and Yudson in Q matrix model IE , and proved later by Wegner for N vector, unitary 
and orthogonal matrix models that a certain class of composite operators with 2s 

fields and gradients endangers the stability of nontrivial fixed points in 2 + e expansions for 
these models. In the following we restrict ourselves to the case of the N vector model. The 
one - loop critical exponents of the corresponding operators in this model read as |2j 

x s = 2s-e S -^^ + 0(e 2 ). (1) 

Judging from these results one could argue that for a sufficiently large s operators become 
relevant (operators with y > 0, y = and y < (y = d — x) are relevant, marginal and 
irrelevant, resp.) and the conventional nontrivial fixed point becomes unstable against an 
infinite number of high-gradient perturbations. This would have important consequences 
for the present understanding of many problems that rely on 2 + e expansions. 

To get a more detailed notion of this problem it seems a reasonable to take advantages 
of the 1/N expansions, which being nonperturbative in its nature relates the 2 + e and 4 — e 
expansions. The latter is commonly believed to be free of any pathologies. Moreover, in this 
approach the corresponding operators look extremely simple — a s , where a is the auxiliary 
field. The critical exponents for this set of operators had been calculated in Ref. ||, and 
independently in Ref. |J: 

I> = 2! + | s(i _ MM^M), (2) 

where 7]\ = 4(2 —ji)Y{2jj, — 2)/T 2 (/i — 1)T(2 — /i)r(/i + 1) and /i = d/2. One can see that the 
problem similar to those in 2 + e expansion arises below d = 3. 

The next attempt to attack the stability problem has been undertaken in the papers , 
where the critical exponents x s have been calculated with the e 2 accuracy. However, the full 
answer for x s given in Refs. J7JQ is not consistent with those (Eq. (fj)) obtained in Refs. []5]|| 
in the 1/N approach, and with the expression for index u (uj = x 2 — d), recently obtained 
in Ref. |TI[ with 1/N 2 accuracy. 

The aim of the present paper is the calculation of the critical exponents x s with 1/iV 2 
accuracy. This provides us the deeper insight into problem and allow to suggest the realistic 
scenario for its solution. The computation scheme developed below is of own interest. We 
remind that only few critical indices are know beyond the 1/N order P-p^]. All of them 



had been calculated by the specific methods, which are not applicable in a general situation. 
With this remark we postpone the discussion to the end of the paper and proceed to the 
calculations. 

The renormalized action for the nonlinear sigma model can be written as : 
S = " y(^) 2 - l*K A M- 2A a + + \oKo. (3) 

Here (p A is the vector field (A = 1, . . . , N), a is the auxiliary scalar field; A and M are 
regularization parameter and mass, respectively. The kernel K is determined from the 
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requirement of the cancellation of the self-energy diagram contribution of the sigma field: 
K{x) = -N/2Gl(x). G^x) is the propagator of the (j) A field (< A (x)0 B (O) >= 5 AB G^x)). 
The regularized kernel is defined as K&(x) = K(x)x~ 2A . The divergencies appearing at 
the calculations of the Feynman diagrams as poles in A are removed by the suitable choice 
of the constants Z\ and Z2. Henceforth we will use the minimal subtraction scheme (MS). 

Unfortunately, the model under consideration is not the multiplicatively renormal- 
ized ||13|| . This means that the freedom in the choice of the renormalization prescriptions 
can not be compensated by the redefinition of the fields and the parameters of Lagrangian. 
Thus one fails to apply the standard RG methods for the calculations of the critical indices. 
The general discussion of this topic can be found in Refs. p~3|,|5f . Here we remind only one 



basic property important for the following, namely, the correlation functions of the fields 
and composite operators in the model under consideration are scale invariant. There is, 
the renormalized one particle irreducible Green's function T° ,n (p, pi) = (0(p)<&(pi)...<&(p n )) 
with operator insertion ($ = <fi, a) satisfies the equation T° ,n (Xp, Xpi) = \ Xo ~ n * Xi >r° ,n (p, P i) . 
As usual, n$x$ means n^x^ + n a x a and the critical dimensions of operator O 

and fields respectively. In the following the composite operators and Green functions are 
assumed to be renormalized, i.e. the all needed counterterms are taking into account. 

The effective algorithm for the computation of the anomalous dimensions of composite 
operators in the first order oil/N expansion exploiting the property of scale invariance has 
been developed in the paper of Vasil'ev and Stepanenko ||. They have shown that the 
simple correlation between the A pole residues of the Green's functions and the correspond- 
ing anomalous dimensions exists. The generalization of this method to all order of 1/N 
expansion (VS scheme) is given below. The basic formula we used for the calculation of the 
anomalous dimensions reads: 

uT°' n (p }Pl ,M) = lim .2Aj2n°T°' n ( P , Pi ,M, A). (4) 

{G} 

Here T°' n (p, pi, M) is the IP I n - point Green's function with operator insertion. A operator 
O is assumed to be the operator with dimension (scaling one). The sum runs over whole set 
of the diagrams (including those with counterterms); is the number of the sigma lines in 
the diagram G; u = —70 + n$7$ (7 is used for anomalous dimensions). 

To derive the Eq. ([|) we note that the scale invariance results in the following form of 
r°' n (p, Pi ,M): 

T > n (p, Pi ,M) = lim J2 T G n (P,Pi,M, A) = p u {M/pfY{p i /p) (5) 

{G} 

where U = 1™ — n<j>x| an . Then acting by MOm on T° ,n (p, Pi,M) and taking into account 
that the only dependence on M in the diagrams results from the propagator of a field 
{G a = M 2A K^ 1 ) (we remind that all counterterms are chosen independent on M (MS 
scheme)) one immediately obtains the Eq. (f|). Note, the finiteness of the lhs of Eq. (|j) 
ensures the cancellation of all A poles in the rhs except for the first order ones. 

Further, in the case of operators {Oi} mixing under renormalization one should seek 
for the proper scaling operators as the linear combinations: Oi = J2 c ikOk- To determine 
both the anomalous dimensions and the form of those one should consider the Eq. (|j) for 
n - point Green's functions with insertion of operators Oi for different n. 
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We apply now the above scheme to the calculation of the critical exponents of a s operators 
in 1/N 2 order. In spite of the fact that there are a lot of operators with the same canonical 
dimension, which could have admixed to a s , this does not happen in the first order of 1/N 
expansion ||. The renormalized a s operator reads as [a s ] = Z s a s (Z s = 1 + q s /NA) and 
is the proper scaling operator in this order. (Henceforth we use the standard notation [O] 
for the renormalized operator.) Going to the next order one finds the counterterms of the 
following form are required - (cr s , a s ~ 2 d 2 cr, d 2 a s ~ 1 ), while all diagrams describing a 2 —>■ <9 4 
transition have not the divergencies. This force us to conclude that in 1/N 2 order the scaling 
operators under consideration have form 

[Os] = [° s ] + ^l* s - 2 d 2 a] + ^PV" 1 ] + O(i). (6) 

To determine ai >2 (really, we need a.\ only) one must consider Eq. (f|) for (s — 1) point 
Green's function of sigma fields with insertion of the operator [O s ]. Doing the same for 
s - point function T 0s,s (p,pi) one obtains the anomalous dimension of the operator in hand. 

However, there is the much simple and elegant way to solve the mixing problem. Let 
us remind that the nonlinear a model is the simplest example of the conformal field theory 
(CFT) (d > 2) P, |14| , |T5|1 . Following along the lines of the paper Jl6| we suppose that the 
conformal operators (CO's) and their total derivatives form a complete basis in the space of 
all operators. In this case any exact scaling operator not being a total derivative of other 
is a conformal one, the opposite is evidently true. (Of course, the more accurate statement 
needs when there is a degeneracy of critical dimensions.) This observation being combined 
with other the well known fact - vanishing of two - point correlator of CO's with different 
scaling dimensions [jHJ - considerably simplifies the solution of the mixing problem. For 
the conciseness we illustrate this idea on the concrete examples, the generalization being 
st raight forward . 

The only nonderivative scaling (and hence conformal) operators on the levels s — 1,2 
(i.e. with the canonical dimensions equal to 2s) have form [0\] = [a] and [0 2 ] = [o" 2 ] + 
a(N)d 2 a. From the requirement of "orthogonality" (02(x)Oi(y)) = one immediately 
obtains ([a 2 (p)]a(-p)) 1PI = a(N)p 2 , witha(iV) = (m/AN) (/i-l)(2/i-3)/(3-/i)+0(l/iV 2 ). 
The above equality holds in all order of 1/N expansion. Note, namely the absence of the 
logarithmic corrections to the above correlator results in the multiplicative renormalization 
of <j s operators in the first 1/N order ||. 

On the level s = 3 two new conformal operators come into a game. On the classical level 
they reads as O3 = a 3 and 0<p = od 2 o — ^"^ c^cr 2 . Then the exact CO's can be written 
as: 

[C ,c,(l)] = [0 (D] + ai[0 W] + + (?a) 

[0 3 c ' (2) ] = [Of] + bxpP] + b 2 d 2 [0 2 ] + hd'lO,] (7b) 

Here a*, bi are some functions of iV and fi. From the "orthogonality" of [O^ 1 ' 2 ^] to [O2], [0\] 
one easy finds that a 2 ~ b 2 ~ 63 ~ 0(1 /N) and 03 ~ 0(l/iV 2 ). In the same time condition 
([03 ] (2?) [03 ] (?/)) = gives the equation entangling bi and a 1; (a x = a®/N + 0(1/N 2 ) 
and 61 = b\ + 0(l/N)): 

([OP](x)[0™ ](y)> + b 1 ([O i l ) }(x)[0?](y)} +a x ([0f ](x)[0f }(y)} = O(N^) (8) 
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To determine the coefficient one must use VS scheme described above (see also Ref. ||): 
6? = 2721/(7^(1) — 7 n (2)); where 7 G i are the anomalous dimensions of the corresponding 

operators in the 1/N order: (All answers it will be given in the units (rji/N) and (r)i/N) 2 
for the first and second order of 1/N expansion, respectively.) 

7 d) = 6(20 - l)[//(2/i - 3) - (/x - l)]/(2 - 0) 

7 n(2 ) = (2/x - l){(2/i - 3) [27^ - 3(110 + 4) + 2 + 90 + 2] - 12(0 - l)}/3(2 - 0). (9) 

The coefficient 721 arises from the mixing of the operators O3 and O3 in 1/N order: 

721 = Jim A£ n<T£' 2 (p,p 1 ,p 2 ) = 20(20 - 3)(40 2 - l)/3(2 - 0) + 0(l/iV 2 ). (10) 
{(?} 

Since 6^ is known, can be easy obtains from the orthogonality condition (Q). One need cal- 
culate two correlators ([0^ 2 ^}(x)[0^' ( ' 2 ^}(y)} in the leading order and one ([O^ ](x)[0^\y)]) 
in the next to leading order. Note, for the determination of the same coefficient in VS scheme 
one need to calculate 20 diagrams in 1/iV 2 order. 

Obviously, the same scheme can be applied for the determination of coefficient ai in the 
operator [O s ] (see Eq. |5|). However, for the persuasiveness we have carried out calculations 
in the both approaches, which result in the same answer a.\ = 70s(s — l)ai: 

(20 - 3)(0 - 1) 4s(0 - 3) + 2 - 50 + 12 

Cti = . 

4(3-0) 4s(20- 3) - 2 - 70+ 12 

Thus this trick allows to fix the form of the scaling operator avoiding the cumbersome 1/iV 2 
order calculations. (The coefficient a% is singular at = S ~ 3/2(1 + l/16(s — 1)). But 
this fact has the simple explanation - at this point the degeneration of 1/N order anomalous 
dimensions occurs (see e.g. formula for b®). So it is only a artifact of used approach, which 
has to be modified in this case.) 

Since the coefficient a\ is known the anomalous dimension of [O s ] can be determined from 
Eq. (fD for T s {p, pi) = i[O s ]{p)a{p\) . . . crip s )). It is instructive to check that the lhs and rhs 
of Eq. (U) have the same momentum dependence. This can be done on the formula level and 
leads us to very nice formula for u^f 1 iuf^ = (— 7^ + 57^)) having obvious resemblance 
with its counterpart in dimensional regularization scheme: 

uf =2A J2n°[KR'G} A . (11) 

G 

Here the sum runs over all diagrams; the KR! operation is the standard operation of the 
subtractions on the divergent subgraphs; [..]a means that only first order poles in A should 
be picked out; is the number of a lines in a diagram G. It is convenient to represent u)P 
as: 

= sis - l)4 2) /2 + sis - - 2) {rf + 2a l7 2i} . (12) 

Here the first two terms arises from the diagrams describing transition a 2 — > a a (iffi), and 
a 3 — > 000 (73 ) in the correlator T^^ip, p^, while the last one is due to the admixture in 
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[O s ]. Since the anomalous dimension of a field (70-) and those of a 2 operator (72 = 2^ a — w 2 ) 
are known up to 1/N 2 order |9| JTT|Jl^] , the problem reduces to the determination of r$ . There 
are 26 diagrams (three 3-loop, ten 4-loop, ten 5-loop, and three 6-loop) which contribute to 
rg . They all can be calculated with the help of technique developed in PJI0|. The final 
expression for 73 reads: 

( 2 ) = /i 2 (/i - l)(15/x 3 - 51/2 2 + 52(1 - 12) 
" 3 = 2(^-2)2 

H(2n - 3) (104// - 383/i 4 + 542/i 3 - 383/i 2 + 104// - 8 



C + 

(13) 



12^-1)^-2)2 

Here C = ip'(l) — ^'(z 1 — 1) and ip(x) = T'(x)/T(x). We do not give the complete expression 
for the critical exponents x s because of its large size. It is more interesting to look on the 
corresponding expansions near 2 and 4 dimensions (we write down only leading in s terms): 



x 2 s +e = 2s- s(se/N) - s/2(se/N) 2 + .. (14a) 
x*~ e = 2s + Qs(se/N) - 3As(se/N) 2 + .. (14b) 

It should be stressed that the coefficient at the last term in Eq. ( |14a| ) differs from those 
obtained in the paper [0. In the 4 — e expansion the exponents x s (corresponding to (4> 2 ) s 
operators) has been so far known with e accuracy only For the additional check of our 
results we have calculated x s up to e 2 order: 

x s = 2s + 66^^ - e 2 s {{s - l)[34(a -2)(N + 8) + 

iv + o 

(lliV 2 + 92iV + 212)] - (13iV + 44)(iV + 2)/2)} /{N + 8) 3 , (15) 

the above expression being in the full agreement with those obtained in the frame of 1/N 
expansion. 

Let us now discuss the obtained results. One can see from the Eqs. fll4]) that as in the 2+e 
as in the 4— e expansions of the exponents x s the second order terms have negative ("wrong" ) 
sign. Thus for a large (se/N) the critical exponents given by the truncated series (|T4|) become 
a negative. This might serve a starting point for the speculations on the stability of the 
conventional fixed points. But, do the first order terms (0) give a well approximation for 
x s when (se/N) is large? The answer, of course, is negative, since the estimative terms 
in (|I4"D are 0((se/N) 3 ) order. Moreover, the pure combinatorical analysis shows that the 
k-th order term in the series ([14]) behaves as s k+1 at large s. Thus to obtain the answer for 
x s which were sensible for large se/N one need sum up all order corrections (see for further 
discussion Refs. fll7|]l8| ). Even taking into account the recent progress in the higher order 



calculations |T2[ the feasibility of this program causes the great doubts. 

The said above concerns both the 2+e, 4— e expansions and the 1/N expansion. However, 
in the case when only a few first terms are calculable, the 1/N expansion has the definite 
advantage in comparison with e ones. Indeed, 1/N expansion is more informative, because 
it contains information on critical exponents in whole interval 2 < fi < 4. In the case under 
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consideration, one can, with some extent of a confidence, judge about the general tendency 
from the first order results. 

Let the function Aj.([i) is the coefficient at the leading in s term (s k+1 ) of the k order term 
in the 1/N expansion of the exponent x s . The function A\(y) and A 2 ([j) are drawn in Fig. [l|. 
It is natural to consider the intervals Ik on which these functions are negative, i.e. if n G Ik 
then Ai < for i = 1, . . . , k. (Indeed, it might be said that the stability problem in the 
2 + e expansion arises from the negativity of the first two coefficients in the expansion ( |14a| ) . 
If one of them were positive, the problem would hardly be considered as serious one.) It is 
seen that the interval I 2 is considerably smaller than I\ ( I\ = [1, 1.5], while I 2 — [1, 1-2]). 
Obviously, if this tendency (Ik —>■ at k —>■ 00) will hold in the higher order, the stability 
problem lost its sharpness. Indeed, in this case for any e > the 1/N corrections starting 
from some order will shift the critical dimensions in the true - irrelevant - direction (at least 
they will be sign - varied). 

We conclude with the following remarks: In all available expansions of the N vector model 
(2 + e, 1/N, 4 — e) there exists the class of operators which acquire a "large" anomalous 
dimensions in the first orders of the corresponding expansions. This fact is not related to 
the expansion used, but has the combinatorical origin. The stability problem arises when 
one tries to extrapolate the first orders results out of their range of applicability. Though we 
do not state the 2 + e expansion is free of the problems [19,20] it seems they are not related 
to the high - gradient operators. 
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FIGURES 




FIG. 1. The functions A\(x) and ^(x) are plotted as the functions of the space dimension. 
The dot line corresponds to A\{x) (the first order in 1/iV), and the solid line — to A2(x) (the 
second order in l/N). 
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